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Abstract 



We derive the next-to-leading order splitting kernels for the scale evolution of fragmenta- 
tion functions for transversely polarized quarks into transversely polarized hadrons. 



Very little is known so far about spin effects in fragmentation of partons into hadrons. An 
exception is semi- inclusive production of A baryons in e + e~, ep, or pp scattering, where the 
polarization of the A has been observed [1-4] and been related to corresponding spin-dependent 
QCD parton-to-A fragmentation functions [5-7]. 

Polarization effects in fragmentation are by themselves interesting since their study opens 
up a new perspective on non-perturbative QCD phenomena in hadron formation. In addition 
it has been realized that, if known with sufficient accuracy, spin-dependent fragmentation func- 
tions can be used as "polarimeters" for nucleon structure [8-11]. For instance, if the polarized 
fragmentation functions for a transversely polarized quark producing a transversely polarized 
A have been extracted accurately in e + e~ annihilation, one can determine with their help the 
much coveted nucleon's transversity densities by studying two spin asymmetries with transverse 
polarization in A production in ep or pp collisions. Behind such a reasoning is of course the 
QCD factorization theorem which states that, for a given produced hadron, the fragmenta- 
tion functions appearing in these scatterings are universal, provided of course the process is 
amenable at all to a description in terms of fragmentation functions. This is the case in the 
presence of a hard scale in the reaction under consideration, such as the virtuality Q of the 
virtual boson in e + e~ annihilation or the transverse momentum of the A in pp — > AX. 

It is expected that future experiments will give very precise information on spin-dependent 
fragmentation functions [12-14]. To analyze such data, an advanced theoretical framework 
is required. According to the factorization theorem the polarized cross section for, say, A 
production in pp^ —> A^X (the arrow denoting transverse polarization) is given by a convolution 
of the form 
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where E,p are the energy and momentum of the produced A. It is again assumed that the 
reaction is characterized by a hard scale such as the transverse momentum pt of the A. We 
denote by s,t,u the Mandelstam variables for the partonic hard-scattering process ab — > cX. 
In Eq. ([[]), f a stands for the unpolarized distribution function of parton a, while the 5fb are 
the transversity densities. The 8D^(z,/j,) represent the transversity fragmentation functions, 
defined in analogy with the transversity distribution functions as 

5D^(z, p) = D$) {Z} n) _ {Zi ^ j (2) 

where D^\z,fi) (D^\z,p,)) denotes the probability for the fragmentation of a transversely 
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polarized parton c to a A with aligned (anti-aligned) transverse spin, carrying the fraction z of 
the parent parton's momentum. Finally, the sum in Eq. (jl|) is over all contributing partonic 
channels a + b — > c + X, with E dba c ah j dp 3 the associated partonic cross section, defined in 
complete analogy with the first line of Eq. ([!]), the transverse polarizations now referring to 
partonic ones: 
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The E dbo c ah j dp 3 are perturbatively calculable thanks to the hard scale involved. 

The factorized form of Eq. (|I|) implies the introduction of a scale fi ~ 0(px), the factoriza- 
tion scale, that reflects the certain amount of arbitrariness in the separation of short-distance 
and long-distance physics embodied in Eq. ([!]). Even though the parton densities and frag- 
mentation functions cannot presently be derived from first principles, their dependence on \i 
is calculable perturbatively in terms of the "DGLAP" evolution equations fl5fl , allowing to 
relate their values at one scale to their values at any other fi. Dependence on /i also arises 
in the procedure of renormalizing the strong coupling constant. Note that in principle one 
could distinguish between factorization scales for the initial and final states and keep also the 
renormalization scale separate; however, for simplicity we keep all scales the same. 

In practice, the precision of the framework of Eq. ([]]) largely depends on the perturbative 
order to which its ingredients are calculated. As already pointed out, there are two places 
where perturbation theory enters: first, the E d8a c ab / dp 3 have the expansion 

«w*i = dMST + (v) di ^ ] + (t)' di ^' + • • • • (4) 

We note that currently partonic cross sections involving transverse polarization in the final 
state are only known at lowest order, except for the reactions e + e~ — > q^q^X and eq" — > eq^X 
for which first-order corrections have been calculated |L6] . 



Secondly, the kernels governing the /x-evolution of the parton densities and fragmentation 
functions also enjoy a perturbative expansion. In this paper we present the first-order correc- 
tions to the evolution of the transversity fragmentation functions defined in Eq. (^|). This seems 
timely in view of forthcoming new experimental information on the production of transversely 
polarized A's [12-14]. In addition, the same evolution kernels [|17| drive the evolution of the 
so-called interference fragmentation functions introduced in Ref. (TO . 



An important observation concerning the 5D^ is that at leading power in the hard scale there 
is no gluonic transversity fragmentation function 5Dg due to angular momentum conservation 
and the helicity-flip nature of the 5D^ [18-20]. This feature also implies that there is no 



2 



mixing with gluons in the /i-evolution of the 5D£, 8D£. The evolution equations are then most 
conveniently written in terms of the combinations 



5D* ± SD} , 



for which they read: 
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Here the superscript "(T)" stands for "time-like" and indicates that we are dealing with a 
fragmentation function. It is instructive to confront Eq. (§) with the corresponding evolution 
equations for the ("space-like") transversity distribution functions 18, [1^, Q in, say, the 
proton, 5q± = 5q± Sq: 
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The evolution kernels 5Pjjq± (U = T, S) occurring in Eqs. (|6|), ([?[) have the perturbative expan- 



sion 
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They are in general not identical, but are closely related to each other, as can be inferred 
from studies of fragmentation functions in the unpolarized case or in the case of longitudinal 
polarization [22-25]. To lowest order, the space-like and time-like splitting functions actually 
do agree, and there is also no distinction between the kernels for the evolution of the + and — 
combinations of densities: 
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where Cp = 4/3, and the -(—prescription is defined in the usual way by 
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The splitting function SPgq± was derived in |26|, . The identity of 8Pgq± and SP^'j}' 

is a manifestation of the so-called Gribov-Lipatov relation (GLR) which connects space- 
like and time-like structure functions within their respective physical regions (£ < 1). The 
space-like and time-like leading order splitting functions are also directly related by analytic 
continuation through x — 1: 
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where x = 1/z and z < 1. Eq. (Jllf) represents the analytic continuation or Drell-Levy-Yan 
relation (ACR) ||29|| . In the second line we have introduced the operation AC that analytically 
continues any space-like function to x — > 1/z > 1 and correctly adjusts the normalization. Note 
that the endpoint contributions oc 5(1 — z) are not subject to the AC operation; however, they 
are necessarily identical in the space-like and time-like cases. 



Both the GLR and the ACR are known to be violated beyond the lowest order |22|, E5 



However, the ACR is based on symmetries of tree diagrams under crossing, and therefore its 
breaking at next-to-leading order (NLO) is essentially of kinematical origin within a given 



regularization prescription, as was shown in Ref. [24]. The main issue here is an extra factor 
of z~ 2e (in MS-scheme dimensional regularization with n = 4 — 2e space-time dimensions) in 
phase space for the time-like case which, when combined with terms singular at e = 0, generates 
extra terms cx \nz in the final answer of the time-like NLO splitting functions. It is then fairly 
straightforward to go through the calculation of Ref. |K| °f the space-like NLO transversity 
splitting functions, and to identify in each contributing Feynman diagram the ACR breaking 
effects in the procedure of analytic continuation to x > 1. 



Combining all extra terms, we obtain in the MS scheme: 
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(for z < 1) where /3 = 11/3 — 4/3T R nf with Ca — 3, T R — 1/2, and the number of active 
flavors, rif. The last term in Eq. ( [12]) obviously represents the breaking of the ACR. It is worth 
pointing out that the structure of that term with its proportionality to both flo and the lowest- 
order splitting function is that of a typical factorization scheme transformation. In other words, 
we could choose a (non-MS) factorization scheme in which the time-like transversity splitting 



functions would be given by AC 
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, without any extra term, so that no breaking 



of the ACR would occur. This possibility was first demonstrated for the unpolarized |3T| and 
longitudinally polarized |24| cases, which are more general in the sense that singlet mixing is 
present there. In the following, we do stay within the more conventional MS scheme, however. 

Inserting the explicit result of 0] for the NLO 5Pqf± (#), performing the analytic con- 



tinuation, and adding the endpoint contributions, we arrive at the final result for 8P qq ± 
We first define 
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where Li2(z) is the dilogarithm function. Then 
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where C(3) « 1.202057. 



We are also in the position now to obtain the difference 5-PiJ± < ' 1 ' ) (-2) — SP^'j} 1 ' \z) which is 
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non-zero if the GLR is violated. Using Eqs. QTB[) and flTT| ) along with the results of Ref. 
we find: 



(z) = C 2 lnztfp^z) (3 + 41n(l - z) -2 In 2) . (18) 



Note that as in the unpolarized and longitudinally polarized cases [T22J, the violation of the 
GLR only occurs in the Cp part of the splitting function. 

For numerical applications it is convenient to have the Mellin-n moments of the NLO split- 
ting functions 5PL ± (z), defined as 
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We obtain: 
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where 77 = ±, and where in the first line we have expressed the result in terms of the moments 
of the space-like NLO transversity splitting functions |30], [32|, The sums appearing above 
are defined by 



n 1 

3=1 J 
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Their analytic continuations to arbitrary Mellin-n (which depend on 77) can be found in |34 

In summary, we have presented the next-to-leading order kernels for the evolution of trans- 
versity fragmentation functions. Our results will become useful in analyses of future precision 
data sensitive to leading-twist transverse-spin effects in fragmentation. 
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